Prom high temperature supercondutivity to quantum spin liquid: progress in strong 

correlation physics 
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This review gives a rather general discussion of high temperature superconductors as an example 
of a strongly correlated material. The argument is made that in view of the many examples of 
unconventional superconductors discovered in the past twenty years, we should no longer be surprised 
that superconductivity emerges as a highly competitive ground state in systems where Coulomb 
repulsion plays a dominant role. The physics of the cuprates is discussed, emphasizing the unusual 
pseudogap phase in the underdoped region. It is argued that the resonating valence bond (RVB) 
picture, as formulated using gauge theory with fermionic and bosonic matter fields, gives an adequate 
physical understanding, even though many details are beyond the powers of current calculational 
tools. The recent discovery of quantum oscillations in a high magnetic field is discussed in this 
context. Meanwhile, the problem of the quantum spin liquid (a spin system with antiferromagnetic 
coupling which refuses to order even at zero temperature) is a somewhat simpler version of the high 
T c problem where significant progress has been made recently. It is understood that the existence 
of matter fields can lead to de- confinement of the 17(1) gauge theory in 2 + 1 dimensions, and novel 
new particles (called fractionalized particles), such as fermionic spinons which carry spin i and no 
charge, and gapless gauge bosons can emerge to create a new critical state at low energies. We even 
have a couple of real materials where such a scenario may be realized experimentally. The article 
ends with answers to questions such as: what limits T c if pairing is driven by an electronic energy 
scale? why is the high T c problem hard? why is there no consensus? and why is the high T c problem 
important? 

PACS numbers: 



I. INTRODUCTION 

High temperature superconductivity was discovered in 
cuprates in 1986. Q The rapid raising of the transition 
temperature to well above the melting point of nitrogen 
shattered the old record of 23 K. Furthermore, the fact 
that high T c superconductivity was discovered in a rather 
unexpected material, a transition metal oxide, made it 
clear that some novel mechanism must be at work. The 
intervening years have seen great strides in high T c re- 
search. The growth and characterization of cuprate sin- 
gle crystals and thin films have advanced to the point 
where sample quality and reproducibility problems which 
plagued the field in the early days are no longer issues. 
At the same time, basically all conceivable experimental 
tools have been applied to the cuprates. Indeed, the need 
for more and more refined data has spurred the develop- 
ment of experimental techniques such as angle resolved 
photoemission spectroscopy (ARPES) and low temper- 
ature scanning tunneling microscopy (STM). Today the 
cuprate is arguably the best studied material outside of 
the semiconductor family and a great deal of facts are 
known. It is also clear that many of the physical proper- 
ties are unusual, particularly in the metallic state above 
the superconductor. Superconductivity is only one as- 
pect of a rich phase diagram which must be understood 
in its totality. 

It is often remarked that there is no consensus for the 
mechanism of high T c superconductivity. This may be 
true but I must emphasize that a lack of consensus is 
not synonymous with a lack of understanding or lack of 



progress. In fact, I will argue that the basic physics of 
the cuprate family is well understood. While there are 
hundreds of high T c compounds, they all share a lay- 
ered structure which contains one or more copper-oxygen 
planes. The low energy physics of these planes can fur- 
ther be simplified to a model of electrons hopping on a 
square lattice called the one band Hubbard model. [1, 3 
(Many of the details left out in this paper can be found 
in a more comprehensive review. 
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where c ia . is the creation operator of an electron with 
spin a on a square lattice, ni a = c\ a Ci a , tij is the hop- 
ping matrix element between sites i and j, (we shall de- 
note nearest neighbor hopping by t and further neighbor 
hopping by t', t" etc.), and U is the repulsive energy cost 
due to screened Coulomb interaction to put two electrons 
with opposite spin on the same site. At half-filling (one 
electron per site), there is a metal-to-insulator transition 
as the ratio U/t is in increased. The insulator is called 
a Mott insulator, Q as opposed to a band insulator, be- 
cause its existence is driven by strong repulsion and is not 
described by band theory. The latter would require the 
state to be metallic at half-filling. The term "strong cor- 
relation physic" is now used to describe the state of affairs 
when interaction energy dominates the kinetic energy in 
controlling the electron motion. As seen in Fig. 1, for 
large enough U/t the electrons prefer to be localized on 
the lattice site because any hopping to reduce the kinetic 
energy t requires double occupation of some site, which 
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FIG. 1: Structure of the Cu-O layer in high T c materials. 
Copper atoms sit on a square lattice with oxygen atoms in 
between. The electronic structure is simplified to a one band 
model shown on the right, with electrons hopping with matrix 
element t. There is an antiferromagnetic exchange J between 
spins on neighboring sites. 



costs U. This insulator is also predicted to be antifer- 
romagnetic (AF), because AF alignment permits virtual 
hopping to gain an energy J = 4t 2 /U by second order 
perturbation theory, whereas hopping is strictly forbid- 
den by Pauli exclusion for parallel spins. It is now gen- 
erally agreed that the parent compound of the high T c 
cuprate is a Mott insulator. 

Things get interesting when electron vacancies (called 
holes) are introduced into the copper-oxygen layers in a 
process called hole doping, i.e. a charge reservoir away 
from the copper-oxygen plane is introduced which re- 
moves electrons for the plane. We denote the concen- 
tration of holes by x. The resulting phase diagram is 
schematically shown in Fig. 2. The AF order is rapidly 
destroyed by a few percent of holes, beyond which su- 
perconductivity appears. The transition temperature 
reaches a maximum around 15% doping, which is called 
optimal doping. The dome-shaped T c is characteristic of 
all hole doped cuprates, even though the maximum T c 
seems to be clustered around two groups. It is about 
40 K in the La^Sr^CuC^ (LSCO) family and 93 K and 
higher in a second family which includes YBa 2 Cu307_,5 
(YBCO) and Ba 2 Sr 2 CaCu 2 9+d - (Bi-2212). The high- 
est T c at ambient pressure of 135 K was reached in 
HgBa 2 Ca 2 Cu30 8+5 in 1993.0 

Before high T c superconductors were discovered in 
1986, almost all superconductors were believed to be 
s-wave BCS superconductors. The requisite attractive 
interaction to form pairs is provided by electrons ex- 
changing a phonon. There were a few potential ex- 
ceptions in a class of strongly correlated metals called 
the heavy fermions, but T c was very low, typically less 
than 1 K.[8j. This is why the discovery of supercon- 
ductivity in a system where the repulsion is strong 
enough to create a Mott insulator was such a surprise. 
However, in the intervening years, we know of many 
examples of non s-wave superconductors. These all 
occur in strongly correlated materials and are clearly 
not driven by electron-phonon coupling. In the heavy 
fermion systems, experimental progress has removed 
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FIG. 2: Schematic phase diagram of high T c materials. The 
antiferromagnet (AF) is rapidly destroyed by doped holes. 
The d-wave superconductor is subject to strong phase fluctua- 
tions below the dotted line, where the proliferation of vortices 
has been detected by the Nernst effect. A pseudogap region 
extends up to high temperatures in the underdoped region. 



any doubt about the non-s-wave nature of a large 
number of compounds. |9( Other new examples include 
Sr 2 RuC>4 exhibiting triplet pairing which breaks time- 
reversal symmetry: [loll hydrated cobaltates,[ll|, and sev- 
eral systems near the quantum critical point of a zero 
temperature ferromagnetic (UGe 2 )or antiferromagnetic 
(CePd 2 Si 2 , Celn3)phase transition. [13, EH The transi- 
tion temperature of these systems remain low, less than 
5 K. Of particular interest to the present discussion is 
the superconductivity discovered in a series of layered or- 
ganic molecular solids. [T^. fl5j| This is because these com- 
pounds also live in the vicinity of the Mott transition. 
Unlike the cuprates, they are not doped and remain at 
half filling. In these materials the effective hopping pa- 
rameter t is sensitive to pressure and to the choice of an- 
ion molecules. The ratio U/t can be tuned right through 
the Mott transition. [l6| Amazingly it was discovered that 
when the Mott insulator is destroyed, the system immedi- 
ately becomes a superconductor, before becoming a metal 
at even higher pressure. Furthermore, the transition tem- 
perature reaches 11.6 K, the highest known among the 
organics. There is also strong evidence that these super- 
conductors have d-wave pairing symmetry. [l7| 

I would argue that 11.6 K for an organic metal qualifies 
it as an example of a high T c superconductor! The rea- 
son is that the electronic energy scale for organic solids 
is much smaller that that for ordinary solids. For ex- 
ample, the hopping matrix element t is about 0.05 eV 
compared with 0.4 eV for the cuprates. Thus the ratio 
ksTc/t w ^ is about the same for both systems. To 
emphasize this point, in Fig 3 I have put both materials 
on the same phase diagram in the parameter space U/t 
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FIG. 3: Location of high T c cuprates and organic supercon- 
ductors in the Hubbard model phase diagram. At half filling, 
the antiferromagnetic insulator onsets when U /t exceeds a 
critical value U c /t, where the Mott transition occurs. High T c 
superconductivity occurs when holes are doped into the Mott 
insulator over a concentration range between 6% and 25%. 
In certain organic compounds, 12 K superconductivity lives 
on the boundary between the Mott insulator and the metal. 
The ratio ksTc/t is about i for both systems. Whether the 
two superconducting regions are connected is not known and 
indicated by the question mark. 



and x of the Hubbard model. Is the d-wave supercon- 
ductor that appears with doping connected with the one 
that appears under pressure? We do not have the answer 
at present. My point is that with so many "unconven- 
tional" examples, our mindset today should be different 
from that of 20 years ago, and we should be more recep- 
tive to the idea that superconductivity may be a highly 
competitive ground state in a pure repulsive model such 
as the Hubbard model. 

With that remark let us return to the cuprates and 
examine the phase diagram in more detail. The region 
between the disappearance of AF and the onset of su- 
perconductivity is complicated by disorder effects, even 
though heroic efforts to make pure samples of YBCO 
have yielded interesting new information. [l8[ We shall 
not discuss this region further. The regions of the phase 
diagram with doping to the left and right of optimal 
is called underdoped and overdoped, respectively. The 
metallic state above T c in the underdoped region has been 
under intense study and exhibits many unusual proper- 
ties not encountered before in any other metal. This is 
shown below the dashed line in Fig. 2 and has been called 
the pseudogap phase. It is not a well-defined phase in 
that a definite finite temperature phase boundary has 
never been found, so the dashed line should be regarded 
as a cross-over. There is now broad agreement that the 



high T c problem is synonymous with that of doping of 
a Mott insulator. It then makes sense to focus on the 
underdoped region, where the battle line between Mott 
insulator and superconductivity is drawn. 

Since we are interested in the case where U is suffi- 
ciently large compared with t for the Hubbard model to 
be in the Mott insulator phase, it is useful to expand in 
t/U. The leading order result is the t-J model 
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The second term is the AF Heisenberg exchange between 
local spins S% = ^c\ a cr a pCip discussed earlier. The non- 
trivial part of the t-J model resides in the projection 
operator P which restricts the Hilbert space to exclude 
the doubly occupied states. The strong Coulomb repul- 
sion now becomes a constraint of no double occupation. 
Compared with the Hubbard model, the Hilbert space 
is reduced from four states per site to three, namely 
spin up, spin down or empty. The parameters of the 
t-J model appropriate for the cuprates is also well es- 
tablished. J ~ 0.13 eV ~ 1500 K, t/J ~ 3 and t'/t is 
negative, of order —0.2, and is believed to vary somewhat 
from compound to compound. [l|| 

Equations (1) and (2) are deceptively simple looking 
Hamiltonians which have defied accurate numerical or 
analytic solution. Nevertheless, the belief among many 
workers in the field is that they contain the rich physics 
of the high T c phase diagram. The situation is not unlike 
QCD, where the Lagrangian is known, but precise un- 
derstanding of confinement and the mass spectrum has 
just begun to emerge from quantum Monte Carlo after 
decades of hard work. To make matters worse, the high 
T c problem at finite doping is analogous to the QCD 
problem with finite quark density, [2(| where accurate nu- 
merical solution is so far not possible due to the fermion 
sign problem. On the other hand, unlike the quark-gluon 
problem the high T c problem has a lot more experimental 
constraint and input. As a result we know a lot about 
the high T c phenomenology which severely limits the the- 
oretical options. 



II. SIMPLE PHYSICAL PICTURE AND THE 
PSEUDOGAP PHENOMENOLOGY 



Let us start with some simple common sense arguments 
to gain some insight into the nature of the problem of a 
doped Mott insulator. Consider a single hole hoping in 
an AF background as shown in Fig. 1. After one hop we 
find a spin surrounded by ferromagnetic neighbors, cost- 
ing an energy of |J from the 3 ferromagnetic bonds if 
the spins are treated as classical S = \. There is a com- 
petition between the exchange energy J and the desire of 
the hole to hop in order to gain the kinetic energy t per 
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hole. For large enough doping the kinetic energy wins 
and we expect a metallic state with some short range AF 
correlation. By comparing xt and J, we expect this to 
onset at x ~ -f ~ ^, in good agreement with the ex- 
perimental finding. This state should be a Fermi liquid 
state. There is a powerful theorem in Landau Fermi liq- 
uid theory commonly called the Luttingcr theorem [2l| 
which states that the area of the Fermi surface is the 
same as that of free fermions, i.e., it is determined by the 
total density of electrons in the unit cell. In our case the 
area is |(1 — x)A BZ where A BZ = (2n/a) 2 is the area 
of the Brillouin zone (BZ). This is exactly what is found 
experimentally. In Fig. 4(d) we show an example of the 
measured Fermi surface. The precise shape can be fitted 
with a hopping model with further neighbor hopping. 

The opposite limit of a few holes (x <C 1) hopping in 
an AF background is less trivial, but by now reasonably 
well understood. The competition with the AF exchange 
causes the effective hopping matrix element to be renor- 
malized downward from t to J.[22,[23|,[24[ The quasipar- 
ticle nevertheless manages to form coherent bands. The 
bands have minima at (± ± j^j.[25j With finite dop- 
ing the Fermi surfaces are ellipses centered at (±5^, ^Ja) 
as shown in Fig. 4(a). Note that the unit cell is dou- 
bled because of AF ordering and the BZ is reduced to 
the diamond in Fig. 4(a). Applying Luttinger theorem 
to the doubled unit cell, the total area of the Fermi 
surface in the reduced BZ is now (1 — x)Arbz where 
Arbz = \Abz- Therefore we conclude that the area of 
each ellipse (hole pocket) is jAbz- Physically it makes 
sense that transport properties are determined only by x 
carriers occupying small Fermi pockets. The theory of a 
few holes in AF background is quite well developed, and 
recently papers applying effective field theory approach 
borrowed from the particle physics literature are partic- 
ularly notable. [5^, [2?| 

We have good understanding of x <C 1 and x > |. 
What happens in between? Here we run into a dilemma. 
We know that AF order is destroyed for x > 0.03, beyond 
which points we have no indication of unit cell doubling. 
If Fermi liquid theory were to hold, what would hap- 
pen to Luttinger theorem? Recall that the nice physical 
picture of small hole pockets rely on the unit cell dou- 
bling. Once that is absent, Luttinger theorem forces us 
to have a "large" F.S., i.e. one with area proportional to 
1 — x. In that case it will be difficult to see how trans- 
port properties will continue to look as if it is given by x 
holes. We note that while the original derivation of the 
Luttingcr theorem was perturbative in the interaction 
strength, the modern derivation by Oshikawa [28l |29| | is a 
topological one and relies on very few assumptions, not 
much beyond the statement that well defined quasipar- 
ticles exist. In principle, the Fermi liquid can develop a 
heavy mass « - so that the conductivity spectral weight 
n/m* x, but experimentally there is no evidence of 
such heavy mass formation. Parenthetically we point 
out that the three dimensional example of doped Mott 
insulator La2- :r Sr :r Ti03 appears to take the heavy mass 
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FIG. 4: (a) Fermi pockets in a doped AF. Dashed lines indi- 
cate the reduced Brillouin zone due to the unit cell doubling 
of the AF. (b) Fermi surface of a tight binding model with 
first and second nearest neighbor hopping, (c) Schematic pic- 
ture of the Fermi arcs, the excitations are gapless when path 
A crosses the arc but are gaped everywhere along path B. 
(d) Experimental data showing the Fermi surface in over- 
doped Tl-2201 (x — 0.25). Colors indicate the intensity of 
low energy excitations. Data from Plate et al. 33j (e) Experi- 
mental data showing the Fermi arc in one quadrant of Figure 
4(c) in underdoped Na2-xCa a ;Cu202Cl2 (x = 0.1). Data from 
K. Shen et al. [34]] 

route. (3fjj 

It turns out that Nature solves this problem in an ex- 
tremely clever and unexpected way. As far as the ground 
state is concerned, the question is moot because it ap- 
pears that once AF is destroyed the system becomes 
superconducting and the Luttinger theorem cannot be 
applied. What about the normal state above the su- 
perconducting T c ? The extensive work using angle re- 
solved photoemission spectroscopy (ARPES) has shown 
that the gapless excitations lie on an arc.[3l|, [32| Any- 
where apart from the arc, the excitations are gapped. 

This situation is sufficiently strange that requires a bit 
more explanation in terms of the experimental observa- 
tion. ARPES measures the spectrum of occupied elec- 
tron states which can be removed by excitations with 
a photon, i.e. it measures the hole spectral function. 
A spectrum is measured at every k point. In a Fermi 
liquid the spectrum consists of a quasiparticle peak at 
energy As one moves along line A in Fig. 4(c), the 
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peak approaches the Fermi energy and disappears as Ek 
crosses the Fermi surface, thus locating its position. This 
is how the Fermi surface in Fig. 4(b) is mapped out. In 
the underdoped case, what happens is that along line 
B, the quasiparticle peak (now quite broad) approaches 
the Fermi energy but stops before crossing it. Instead it 
loses weight and disappears. You might say that this also 
happens in Fig. 4 (a), if path B misses the hole pocket. 
The important difference is that in the case of the hole 
pocket where there is unit cell doubling, in the extended 
zone scheme we expect the occupied band to live also out- 
side the reduced BZ. Thus if we follow path A we should 
see a quasiparticle peak appearing at the second crossing 
of the ellipse. Further along path A this peak will then 
move down in energy away from the Fermi energy. In 
other words, the back side of the hole pocket should be 
visible in ARPES as an occupied quasiparticle state rises 
up to meet the FS. The surprise is that in Fig. 4(c), there 
is no back side to the ellipse, and one is left with what is 
called a Fermi arc. The spectrum is gapped everywhere 
except on the arc, and the gap reaches a maximum near 
(71", 0). In the superconducting state the arc also becomes 
gapped, leaving a single gapless point along (w, tt) called 
the nodal point. In this way the quasiparticle spectrum 
of a d wave superconductor is smoothly formed out of the 
pscudogap normal state. The size of the maximum gap 
has been measured as a function of doping, and is found 
to increase with decreasing x, in a way which tracks the 
onset of the pseudogap phase shown in Fig. 2. (A word of 
caution: the electron spectral function near (n, 0), the so 
called anti-nodal direction, is extremely broad and, un- 
like the nodal directions, does not show a quasiparticle 
peak even in the superconducting state in strongly under- 
doped samples. So the gap is measured as the pull-back 
of the leading edge of the spectrum for the Fermi en- 
ergy. This picture is corroborated by extensive scanning 
tunneling microscopy (STM) work, which has excellent 
energy and spatial resolution, but no momentum space 
information. 

The Fermi arc and anti-nodal gap is one part of 
a rich phenomenology associated with the pseudogap 
phase. The first evidence of an energy gap comes from 
NMR data, which found that in underdoped samples the 
Knight shift, which measures the spin susceptibility, does 
not follow Pauli's prediction of being temperature inde- 
pendent, but starts dropping around room temperature 
and has lost about 80% of its value by the time T c is 
reached. [1^] The gap also shows up in c-axis frequency 
dependent conductivity, [36| where electrons move from 
one layer to the next, but the transport within the plane 
remains metallic, being dominated by the quasiparticles 
near the Fermi arc. 

What is the origin of the large gap at (71", 0)? One 
suggestion is that it is simply a d-wave superconducting 
gap and the pseudogap phase should be understood as a 
superconductor destroyed by strong phase fluctuations. 
Phase fluctuation is controlled by the superfluid stiffness 
which is given by the superfluid density. Thus T c should 



be proportional to the superfluid density, i.e. decreases 
with decreasing x.JStJ This picture is made more quan- 
titatively accurate if the reduction of the superfluid den- 
sity due to thermal excitations of nodal quasiparticles is 
taken into account. [51, |38| The opposite trend of the en- 
ergy gap and T c is explained, but the origin of such a 
large superconducting gap close to an insulator becomes 
even a greater puzzle. Experimentally this is now ample 
evidence that the picture of fluctuating superconductiv- 
ity indeed survives up to much higher temperature com- 
pared with a conventional metal: perhaps 2 or 3 times 
T c . However, that scale decreases with decreasing dop- 
ing and does not reach up to the pseudogap temperature. 
The key experiment that mapped out this region is the 
Nernst effect, which is sensitive to mobile superconduct- 
ing vortices. [39( For this reason we refer to this region as 
the Nernst region in Fig. 2. 

If the pseudogap is not a pairing gap, it presents a 
great challenge for theory, because while the Fermi arc 
scenario interpolates between the small hole pocket and 
large Fermi surface beautifully, it is not allowed by con- 
ventional band theory or Fermi liquid theory. Fermi sur- 
faces do not simply terminate. This is a situation not 
encountered before in solid state physics. This is why 
the psuedogap phenomenon is considered one of the cen- 
tral mysteries of the high T c story. We shall return to 
this issue in section V. 



III. THE RVB PICTURE 

While the pseudgap phenomena is really strange, to 
a large extent it was anticipated by theory. Here I am 
referring to the concept of the resonating valence bond 
(RVB) introduced by P.W. Anderson[3j and the slave bo- 
son mean field theory and elaborations which followed. 
Here we provide a brief review. 

We explained in the last section that the Neel spin or- 
der is incompatible with hole hopping. The question is 
whether there is another arrangement of the spin which 
achieves a better compromise between exchange energy 
and the kinetic energy of the hole. For S — h it ap- 
pears possible to take advantage of the special stability 
of the singlet state. The ground state of two spins S 
coupled with antiferromagnetic Heisenberg exchange is a 
spin singlet with energy —S(S+ 1) J. Compared with the 
classical large spin limit, we see that quantum mechan- 
ics provides an additional stability in the term unity in 
(S + 1) and this contribution is strongest for S = \. Let 
us consider a one-dimensional spin chain. A Neel ground 
state with S z = ±i gives an energy of — \J per site. 
On the other hand, a simple trial wavefunction of singlet 
dimers already gives a lower energy of — | J per site. This 
trial wavefunction breaks translational symmetry and the 
exact ground state can be considered to be a linear super- 
position of singlet pairs which are not limited to nearest 
neighbors, resulting in a ground state energy of 0.443 J. 
In a square and cubic lattice the Neel energy is — ^ J and 
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— | J per site, respectively, while the dimer variational en- 
ergy stays at — | J. It is clear that in a 3D cubic lattice, 
the Neel state is a far superior starting point, and in two 
dimensions the singlet state may present a serious com- 
petition. Historically, the notion of a linear superposition 
of spin singlet pairs spanning different ranges, called the 
resonating valence bond (RVB), was introduced by [4(| 
and [41| as a possible ground state for the S = \ antifer- 
romagnetic Heisenberg model on a triangular lattice. The 
triangular lattice is of special interest because an Ising- 
like ordering of the spins is frustrated. Subsequently, it 
was decided that the ground state forms a \/3 x y3 su- 
perlattice where the moments lie on the same plane and 
form 120° angles between neighboring sites. [42jl 

Soon after the discovery of high T c superconductors, 
Anderson Q revived the RVB idea and proposed that 
with the introduction of holes the Neel state is destroyed 
and the spins form a superposition of singlets. The va- 
cancy can hop in the background of what he envisioned 
as a liquid of singlets and a better compromise between 
the hole kinetic energy and the spin exchange energy 
may be achieved. Many elaborations of this idea fol- 
lowed, but here we argue that the basic physical picture 
described above gives a simple account of the pseudo- 
gap phenomenon. The singlet formation explains the de- 
crease of the uniform spin susceptibility. The vacancies 
are responsible for transport in the plane. The conduc- 
tivity spectral weight in the ab plane is given by the hole 
concentration x and is unaffected by the singlet forma- 
tion. On the other hand, for c-axis conductivity, an elec- 
tron is transported between planes. Since an electron 
carries spin i, it is necessary to break a singlet. This 
explains the gap formation in <j c (uj) and the energy scale 
of this gap should be correlated with that of the uni- 
form susceptibility. In photoemission, an electron leaves 
the solid and reaches the detector, the pull back of the 
leading edge simply reflects the energy cost to break a 
singlet. 

A second concept associated with the RVB idea is the 
notion of spinons and holons, and spin charge separa- 
tions. Anderson postulated that the spin excitations in 
an RVB state are S — \ fermions which he called spinons. 
This is in contrast with excitations in a Neel state which 
are S — 1 magnons or S — gapped singlet excitations. 

Initially the spinons are suggested to form a Fermi 
surface, with Fermi volume equal to that of 1 — x 
fermions. (43l| Later it was proposed that the Fermi sur- 
face is gapped to form d-wave type structure, with maxi- 
mum gap near (0, 7rh[44| This k dependence of the energy 
gap is needed to explain the momentum dependence ob- 
served in photoemission. 

The concept of spinons is a familiar one in one- 
dimensional spin chains where they are well understood 
to be domain walls. In two dimensions the concept is a 
novel one which does not involve domain walls. Instead, 
a rough physical picture is as follows. If we assume a 
background of short range singlet bonds, forming the so- 
called short-range RVB state, a cartoon of the spinon is 





(a) 



(b) 



FIG. 5: A cartoon representation of the RVB liquid or sin- 
glets. Solid bond represents a spin singlet configuration and 
circle represents a vacancy. In (b) an electron is removed from 
the plane in photoemission or c-axis conductivity experiment. 
This necessitates the breaking of a singlet. 



shown in Fig. \5\ If the singlet bonds are "liquid," two 
S = 5 formed by breaking a single bond can drift apart, 
with the liquid of singlet bonds filling in the space be- 
tween them. They behave as free particles and are called 
spinons. The concept of holons follows naturally [45| as 
the vacancy left over by removing a spinon. A holon 
carries charge e but no spin. 



IV. PROJECTED WAVEFUNCTION, SLAVE 
BOSON AND THE GAUGE THEORY 
FORMULATION OF THE RVB PICTURE 

Is there any calculated tool or mathematical formal- 
ism to put some meat into the physical picture of RVB 
described in the last section? As far as computation is 
concerned, the use of the projected wavefunction has en- 
joyed considerable success. The idea is to write down a 
trial wavefunction of the type 



* = Pr, 



(3) 



where Pq = 1^(1 — n %\ n i\) is called the Gutzwiller pro- 
jection and <j) is any Hartree Fock or BCS wavefunc- 
tion, usually suggested by mean field theory described 
below. The role of the Gutzwiller projection is to re- 
move all doubly occupied states in <f>. Equation (3) is a 
suitable variational wavefunction for the t-J model be- 
cause the constraint is satisfied by definition, and its ex- 
pectation values and correlation functions can be com- 
puted by efficient Monte Carlo algorithms. [irj [47[ The 
mean field parameters can be treated as variational pa- 
rameters. The projection wavefunction gives excellent 
ground state energy and sublattice magnetization at half- 
filling, capturing the important quantum fluctuations of 
the Neel ordered state. With doping it predicted cor- 
rectly the d-wave pairing ground state. [4a| even though 
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the prediction of the co-existence of superconductivity 
with AF up to x rs 0.11 is not in agreement with ex- 
periment. Putting aside the question of whether it is 
the ground state, a comparison of the physical proper- 
ties of the projected d-wave pairing states with a variety 
of experiments was successfully made.Ji^] The trial wave- 
function can be further improved by the Lanczos method 
of repeatedly hitting it with the Hamiltonian. There is 
some controversy as to whether the ground state of the t- 
J model with nearest-neighbor hopping only is a d-wave 
superconductor. (Hoi l5l| but it is clear that superconduc- 
tivity is a highly competitive state, as found by other 
numerical methods such as density matrix renormaliza- 
tion group, cluster dynamical mean field theory [53| 
and variational cluster approximation. [54l [55[ Recently 
it was found that introduction of t' considerably stabi- 
lizes the d-wave superconducting state. [56[ At present, I 
would stay that there is strong numerical evidence that 
the d-wave superconductor is a strong contender for the 
ground state of the t-J model. 

What about analytic theory and where does the mean 
field 4> come from? A useful method is called the 
Gutzwiller approxmation, which imposes the constraint 
approximately by treating the available confi gur ation for 
hopping and exchange in a statistical basis. [53] This is 
clearly related to the slave-boson method which we dis- 
cus below. The slave-boson method was developed for 
the Kondo problcm[58, 59] It has enjoyed great success as 
the best way to undersand the properties of a remarkable 
class of materials called the heavy fermion compounds, 
where Fermi liquid theory has been stretched to the ex- 
treme, with effective mass as large as several thousand 
times the free electron mass. 60| The idea is to write the 
electron operator as a product of the boson and fermion 
which carries the spin index 



with the condition 



1. 



(4) 



(5) 



This constraint can be enforced with a Lagrangian mul- 
tiplier Ai. Note that Eq. (4) is not an operator identity 
and the right-hand-side does not satisfy the fermion com- 
mutation relation. Rather, the requirement is that both 
sides have the correct matrix elements in the reduced 
Hilbert space with no doubly occupied states. For exam- 
ple, the Heisenberg exchange term is written in terms of 
fL ha only @ 
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We then decouple the exchange term in both the particle- 
hole and particle-particle channels via the Hubbard- 
Stratonovich (HS) transformation. 

Then the partition function is written in the form 



Z= / DfDpDbDXDxDA exp 



drLi 



(7) 
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with Xij representing fermion hopping and Aij repre- 
senting fermion pairing corresponding to the two ways 
of representing the exchange interaction in terms of the 
fermion operators. J = 3J/8 is chosen to reproduce the 
mean field self-consistent equation which is obtained by 
the Feynman variational principle. 

Mean field theory corresponds to the saddle point solu- 
tion to the functional integral. The mean field conditions 
are 



Xij — J2a(ficrfjv) 
&ij = (filfjl - filfit) 



(6) 



(9) 
(10) 

Let us write Xij = \Xij\ e%aii an d ignore the fluctuation 
of the amplitude. Furthermore, in the last term in Eq. (8) 
we replace fj a fja by Xij- Then the rather complicated 
Lagrangian of Eq. (8) has a rather simple interpretation. 
It describes fermions and bosons hopping on a lattice 
with hopping matrix element Xij ■ The phase lives on 
the links, and plays the role of the spatial components of 
a lattice gauge field, while the fields introduced to en- 
force the constraint become the time components. Note 
that both fermions and bosons are coupled to the same 
gauge field. In addition, the fermions may have singlet 
pairing amplitude given by Ay. Thus we come to the 
conclusion that the t-J model is equivalent to a lattice 
gauge theory with non-relativistic fermions and bosons 
coupled with a compact U(l) gauge field (compactness 
simply refers to the fact that the gauge field ay is a 
phase defined module 2tt). At this point the gauge field 
has no dynamics. It fluctuates freely and is in the infinite 
coupling limit. Interesting dynamics emerge upon inte- 
grating out some of the matter field but one is left with a 
lattice gauge theory with strong coupling. The mapping 
is basically exact, but the question remains as to how to 
deal with such a model. 

Before discussing the importance of gauge fluctuations, 
let us examine some examples of mean field solutions. 
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FIG. 6: Schematic phase diagram of the Z7(l) mean field 
theory. The solid line denotes the onset of the uniform RVB 
state (x / 0). The dashed line denotes the onset of fermion 
pairing (A 7^ 0) and the dotted line denotes mean field Bose 
condensation (6 / 0). The four regions are (I) Fermi liquid 
X 0, b 0; (II) spin gap X / 0, A / 0; (III) d-wave 
superconductor x 7^ 0, A / 0, b ^ 0; and (IV) strange metal 
X / 0. From Lee and Nagaosa.[6ll] 

1. d-wave pairing states. 

Here Xij — X 1S constant and A™ = A for (ij) 
bonds along x and Ay = —A for (ij) along y, i.e. 
it has d-wave pairing geometry. Without pairing, 
the fermions hop on a tight binding band with dis- 
persion 

£/(fc) = — 2 Jx(cos k x a + cos k y a). (11) 
With pairing we have the classic <i-wave dispersion 

E(k) = ^f(e f (k)-n f ) 2 + \A k \i (12) 

where [if is the fermion chemical potential and 
At = 2A(cosk x a — cosk y a). The bosons see the 
same band dispersion and condense at the bot- 
tom of the band minimum at low temperatures. 
In mean field theory the Bose condensation tem- 
perature is proportional to the boson density x. 
Below this temperature we have electron pairing, 
because the BCS order parameter (c^c^ki) — 
^o(/fct/-fel) where (b) = b Q . The mean field 
phase diagram is shown in Fig. 6 and captures some 
key features of the high T c phase diagram shown in 
Fig. 2. In particular, the d-wave superconducting 
state appears at intermediate doping and a spin 
gap state (region II) where a d- wave-like gap exists 
for spin excitations but not charge excitations, an- 
ticipates many of the properties of the pseudogap 
phase. 

2. The staggered flux state. 

Early in the development of the mean field theory, a 




FIG. 7: The energy dispersion of the staggered flux 
phase. Note the massless Dirac spectrum at the nodal points 
(±5, ±5). Figure shown is for the special case of 7r-flux. In 
general, the nodal spectra becomes anisotropic. With doping 
Fermi pockets are formed when the Fermi energy crosses the 
energy spectrum. 



variety of mean field states were discovered which 
give identical dispersion. Notable among these is 
the staggered flux state [13] In this state the hop- 
ping Xij is complex, Xij = Xo exp (i(-l) l -+^$ ), 
and the phase is arranged in such a way that it 
describes free fermion hopping on a lattice with a 
fictitious flux ±4$o threading alternative plaque- 
ttes. Remarkably, the eigenvalues of this problem 
are identical to that of the ci-wave superconductor 
given by eq. (12), with fif — and 

tan$ = — . (13) 

X 

The case <f>o = 7r /4, called the ir flux phase, is 
special in that it does not break the lattice trans- 
lation symmetry. As we can see from Eq. (13), 
the corresponding <i-wave problem has a very large 
energy gap and its dispersion is shown in Fig. 7. 
The key feature is that the energy gap vanishes at 
the nodal points located at (±f , ±f ) ■ Around the 
nodal points the dispersion rises linearly, forming 
a cone which resembles the massless Dirac spec- 
trum. For the ir flux state the dispersion around 
the node is isotropic. For $0 less than 7r/4 the gap 
is smaller and the Dirac cone becomes progressively 
anisotropic. The anisotropy can be characterized 
by two velocities, vp in the direction towards (ir, ir) 
and da in the direction towards the maximum gap 
at (0,tt). 

The reason various mean-field theories have the 
same energy at half-filling was explained by Affleck 
et a/. [63[ and Dagotto et aZ.[64| as being due to a 
certain SU (2) symmetry. It corresponds to the fol- 
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lowing particle-hole transformation 

fh - mfh+Pifu (14) 

Note that the spin quantum number is conserved- 
It describes the physical idea that adding a spin-up 
fermion or removing a spin-down fermion are the 
same state after projection to the subspace of singly 
occupied fermions. It is then not a surprise to learn 
that the Gutzwillcr projection of the d-wave super- 
conductor and that of the staggered flux state gives 
the same trial wavefunction, up to a trivial overall 
phase factor, provided fj,f — and Eq. (13) is sat- 
isfied. A simple proof of this is given by Zhang et 
ai.[13|. The energy of this state is quite good. The 
best estimate for the ground state energy of the 
square lattice Heisenberg antiferromagnet which is 
a Neel ordered state is (Sj -Sj) = —0.3346 J. [65L l66l| 
The projected it flux state [48| gives — 0.319J, which 
is excellent considering that there is no variational 
parameter. 

In the presence of a hole, the dispersion of the 
staggered-flux phase is still given by Eq. (12) with 
fXf = 0, but the Fermi level is now located at fif < 
which lies below the node. This is shown in Fig. 7. 
If the bosons are condensed, this becomes a Fermi 
liquid state with small hole pockets just like what 
is shown in Fig. 4(a). This state has higher energy 
than the d-wave superconductor because in a su- 
perconductor the node is shifted in k space away 
from (n/2a, n/2a) towards the origin, but its en- 
ergy is always tied to the Fermi level. The stag- 
gered flux state was proposed by Hsu, Marston and 
Affleck[67] in 1991 to be the origin of the pseudo- 
gap state. They pointed out that with doping, the 
staggered flux state has the remarkable property 
that orbital currents flow around each square pla- 
quette in a staggered way, so that the physical or- 
der parameter of this state is the staggered orbital 
current. This proposal did not receive much at- 
tention because the appearance of an ordered state 
with Ising order parameter requires a finite tem- 
perature phase transition, which has never been 
seen experimentally. Furthermore, the model re- 
quires hole pockets instead of the arcs which were 
later observed. While the matrix element effect 
can reduce the spectral weight of the backside. [68j 
the model requires sharp quasiparticle peaks which 
should be observable, especially near the end of the 
arc, where the weight reduction is only a factor of 
2. In 2002, Chakravarty et al.\&^ revived this pro- 
posal, arguing that disorder effect may round the 
transition. They named this state the d-density 
wave (DDW) state. Since this state is in fact iden- 
tical to the staggered flux state of Hsu et al.[67], 
the two names are sometimes used interchangeably 
in the literature. However, the philosophy of their 




doping (x) 

FIG. 8: Energy per site of projected wavefunctions in units 
J. SF stands for staggered flux. From Ivanov. [74| 

approach is very different. Chakravarty et al. take 
a phenomcnological Landau theory approach. For 
them the superconducting state is an entirely dif- 
ferent order parameter and the two orders compete 
with each other. I believe this view misses the key 
fact that all this is happening close to the Mott 
insulator. While the d-wave superconductor and 
the staggered flux (DDW) states are very different 
states at the mean field level, the mean field states 
do not respect the no double occupation constraint. 
The story changes if we enforce the constraint by 
applying the Gutzwillcr projection. As mentioned 
earlier, the projected states become identical at half 
filling and by continuity must share a lot of similar- 
ity slightly away from half filling. For example, we 
found that at x = 0.1 the projected d-wave super- 
conductor has short-range orbital current order pTOj 
and presumably the projected staggered flux state 
has short range d-wave superconductor order. Both 
have significant short range AF order [ti| • This 
kind of consideration motivated Wen and I to in- 
troduce the SU (2) gauge theory in 1996 as an im- 
provement over the Hsu et al. proposal. [72} Instead 
of an ordered state, the pseudogap phase is consid- 
ered to be a fluctuating state which included the 
staggered flux state and d-wave superconductivity 
on equal footing. This will be described in more 
detail in the next section. 

It is interesting to examine the energetics of various 
projected states as shown in Fig. 8. It was found 
that the best state is a projected d-wave supercon- 
ductor and the sublattice magnetization is nonzero 
for x < x c , where x c = 0.11 for t/J = 3. [74] The 
projected staggered-flux state always lies above the 
projected d-wave superconductor, but the energy 
difference is small and vanishes as x goes to zero, 
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as expected. The staggered- flux state also prefers 
antifcrromagnetic order for small x, and the critical 
F is now 0.08, less than that for the projected d- 
wave superconductor. The projected staggered flux 
state is the lowest energy nonsuperconducting state 
that has been constructed so far. 



SU(2), VORTEX CORE AND THEORY OF 
THE PSEUDOGAP PHASE 



The lesson from examining projected wavefunctions is 
that the close relationship between different mean field 
states must be taken seriously. In fact, the d-w&ve state 
and the staggered flux state are not the only states which 
have similar energies after projection. There is a con- 
tinuous family of states related by the SU (2) symmetry 
given by Eq. (14). Wen and I developed a formulation to 
take this into account by introducing an SU(2) doublet 
of bosons (&i, 62) instead of the single boson in the 17(1) 
gauge theory. [73] For our purpose this is just a technical 
way of generating different mean field states which are 
parametrized by rotating the quantization axis I in the 
SU(2) space. [73 This way different states can be visual- 
ized and smoothly connected to each other in space and 
time. In Fig. 9 we show such a representation. The north 
and south poles correspond to the two degenerate stag- 
gered flux (DDW) states which break translational sym- 
metry. The equator is the c£-wave superconducting state. 
In between are many states which share both kinds of or- 
der. One advantage of this approach is that we can con- 
struct a model of the vortex core of an hc/2e vortex. (75j 
This takes the form of a meron, or half a skyrmion. As 
shown in Fig. 10 the center of the vortex is occupied by 
the staggered flux state. This solves a serious problem 
with the original £7(1) formulation which favors the hc/e 
vortex [7a, ul\ because it is energetically favorable to re- 
tain the fermion pairing and make the boson wind by 
27r around the vortex. That approach ignored the fact 
that another state, the staggered flux state, is available 
to take the place of the normal state inside the core, al- 
lowing us to construct an hc/2e vortex that costs very 
little energy. 

Our model of the vortex core also explains a puz- 
zling experimental observation, i.e. STM tunneling found 
that an energy gap remains when tunneling into the core 
region. 78] This is opposite to what was found in a con- 
ventional superconductor, where bound states are formed 
inside the core which fills in the gap. j79j 

The SU(2) construction allows us to construct a sim- 
ple picture of the finite temperature phase diagram 
which accounts for the Nernst region and the pseudogap 
phase. [75[ This is shown in Fig. (11). At low tempera- 
ture the SU{2) quantization vector J is in the x-y plane 
and we have the d-wave superconductor ground state. As 
the temperature is raised, vortices are created with stag- 
gered flux cores and the proliferation of these vortices 
drives the Berezinskii-Kosterlitz-Thouless transition in 
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FIG. 9: The quantization axis J in the 5(7(2) gauge theory. 
The north and south poles correspond to the staggered flux 
phases with shifted orbital current patterns. All points on the 
equators are equivalent and correspond to the d-wave super- 
conductor. 



the standard way. Above this transition (which is the 
true superconducting T c ), vortex and anti- vortex prolif- 
erates giving rise to the Nernst phase. At even higher 
temperature the I vector is completely disordered and 
this is our picture of the pseudogap phase. 

We emphasized that in this theory the pseudogap 
phenomenon and superconductivity are intimately con- 
nected. There is no separate pairing mechanism for su- 
perconductivity. What drives superconductivity is the 
coherence of the boson to select the true ground state 
out of a myriad of fluctuating possibilities. This is very 
different from the competing order scenario which does 
require a separate pairing mechanism and a completely 
separate energy gap scale. This dichotomy has spurred a 
debate concerning one gap vs. two gaps, i.e. whether a 
smaller energy gap appears which scales with T c .[8(| The 
debate is often couched in a black and white language, 
with one gap favoring a superconducting gap destroyed 
by phase fluctuations, and two gaps implying the need 
for some kind of competing order. In my view, the truth 
is likely to be more complicated. In the mean field RVB 
picture, a sharp quasiparticle peak appears below T c with 
weight x which follows the d-wave dispersion with only 
one gap. The mean field picture is probably too sim- 
plistic. For example, it is possible that the large gap at 
(0, 7r) as a spin gap which can remain broad while the 
low energy quasiparticles near the nodes become coher- 
ent below T c . To a first approximation, the coherent 
nodal quasiparticle has a dispersion which extrapolates 
to the large pseudogap at (0,7r). There could be a co- 
herence energy scale which scales with T c , but exactly 
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FIG. 10: Model for a hc/2e vortex. The SU(2) quantiza- 
tion axis points for the north pole at the center (forming a 
staggered flux vortex core) and rotates smoothly towards the 
equatorial plane as one moves out radially. 



how this coherence scale affects the density of states and 
how it develops as a function of temperature is an open 
question. 

The issue of one gap/two gaps is not settled experimen- 
tally. For moderately underdoped Bi2212 (T c > 50 K) 
the evidence from ARPES is that the quasiparticle in 
the superconducting state is reasonably peaked even near 
the antinodes and obey the d-wave dispersion with a sin- 
gle gap which increases with decreasing doping. This 
is supported by low temperature thermal conductivity 
data which measures the ratio Vf/va where is the 
quasiparticle velocity in the direction of (0, 7r) . [8l|] It is 
found that v& increases with decreasing x and extrap- 
olates to the antinodal gap measured by ARPES. On 
the other hand, for severely underdoped samples and one 
layer cuprates with low T c 's, there are claims based on 
ARPES that the energy gap in the Fermi arc region near 
the nodal point does not scale with the pseudogap at 
(0,7r), which increases with decreasing doping as men- 
tioned before. [82| Instead, it seems to stay constant or 
increase with decreasing doping. It is argued that this 
reveals a new gap scale associated with superconductiv- 
ity. I should caution that deeply underdoped samples 
are known to be strongly disordered, and the disorder in- 
creases with reduced doping. Furthermore, the lineshape 
remains very broad in the antinodal direction even in the 
superconducting state. Thus it is risky to draw strong 
conclusion from lineshapes without an understanding of 
disorder effects and of the lineshape. 

Other support for two gaps comes from Andreev 
reflection studies [S3l] and Raman scattering. (84] In a 
superconductor-normal metal junction in conventional 
superconductors, normal electrons incident on the junc- 
tion has an extra channel for transport, by tunnelling 
as a Cooper pair into the superconductor and Andreev 
reflected as a hole. This leads additional conductance 
below an energy scale of the energy gap. Such extra con- 
ductance was observed in underdoped cuprates, but the 
energy scale observed is much lower than the pseudogap 
and is more related to T c . I note that in contrast to con- 
ventional tunnelling, Andreev reflection does not simply 
measure the density of states, but requires coherence of 
the quasiparticle in its interaction with the condensate. 
What is seen in the Andreev data may be this coherence 
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FIG. 11: Schematic picture of the quantization axis I in dif- 
ferent parts of the phase diagram shown in Fig. 2. (a) In the 
superconducting phase J is ordered in the x-y plane, (b) In 
the Nernst phase, J points to the north or south pole inside 
the vortex core, (c) The pseudogap corresponds to a com- 
pletely disordered arrangement of i". (i" is a three dimensional 
vector and only a two dimensional projection is shown.) 



scale. 

I must emphasize that the simple cartoon shown in 
Fig. 11 is only an approximate picture. We have assumed 
that the bosons are locally condensed and can be treated 
as a c- number which varies in space and time. However, 
even at T = 0, the vortex configurations shown in Fig. 10 
can tunnel between each other and destroy the staggered 
flux order at some time scale. I think the correct answer 
requires a quantum mechanical treatment of the boson 
strongly coupled to gauge fields, which is not available 
at present. In particular, we have not yet been able to 
compute the ARPES spectrum and make a satisfactory 
comparison with experiment. We make crude approxi- 
mations such as assuming a binding of the bosons with 
fermions via gauge fluctuations. (73j As an example of an 
alternative aproach, Ribeiro and Wen [85[ introduced a 
new formulation which hybridizes the physical hole with 
the spin carrying fermions and have had success in under- 
standing the higher energy spectra. Their theory seems 
to favor the two gap scenario. The truth is that the 
theory of the spectral function is not under control at 
present: the problem of a fermion and boson strongly in- 
teracting with a gauge field at finite temperatures is too 
daunting with currently known techniques. Rather than 
trying to explain data in detail, it may be more fruitful 
to step back and attempt a classification of the pseudo- 
gap state. We have developed a view that the pseudogap 
phase belongs to the deconfined side of the lattice gauge 
theory, so that discussion in terms of fermions, bosons 
and noncompact gauge fields makes sense. More explic- 
itly, we propose that the pseudogap is best understood as 
the doping of a particular spin liquid called the algebraic 
spin liquid. This is explained in the next section. 
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VI. GAUGE THEORY, DE-CONFINEMENT 
AND SPIN LIQUIDS 

Let us return to the gauge theory formulation de- 
scribed in section III. As it stood Eq. (8) is pretty much 
an exact reformulation of the t-J model. The constraint 
is enforced exactly upon integrating out the gauge fields. 
However, the gauge field is fluctuating strongly because 
there is no restoring force, and critics of the gauge theory 
approach have pointed to this as evidence that this ap- 
proach cannot be trusted to give meaningful results. At 
a deeper level, I think what lies behind this objection is 
the implicit assumption that a phenomenon called con- 
finement, familiar in QCD, necessarily takes place when 
the gauge coupling is large. In this case, fermions are 
tightly bound to bosons by the gauge field, and one just 
recovers the original intractable electron problem with 
constraint. In contrast, our way of thinking assumes that 
confinement does not happen, in which case the fermions, 
bosons and gauge fields retain their identity, albeit with 
strong interaction. The ultimate strongly coupled fixed 
point is nontrivial and different from that of free parti- 
cles, but as is usual in the strong coupling problem, it 
may be possible to access their behavior using artificial 
expansion parameters such as j* where N is the number 
of copies of the matter field. Even though the physical 
problem may correspond to N ~ 2 or 4, as long as con- 
finement does not take place the behavior of the physical 
system is smoothly connected to the large N system and 
the physical behavior of the system can be understood, 
even though the critical exponents cannot be computed 
quantitatively. 

In the case of a compact U(l) gauge, it is well known 
that in the 2+1 dimension the theory is confining even 
for arbitrarily small coupling, due to the proliferation of 
instantons. The point is that in a compact gauge field, 
a 27r flux can pop up through an elementary plaquette. 
The space-time point where this happens is called an in- 
stanton event, which is also called a magnetic monopolc 
in Euclidean space-time. If these monopoles proliferate, 
the total flux through the sample is not conserved. The 
strongly fluctuating gauge flux leads to confinement be- 
tween external charges. At first sight, this bodes poorly 
for the gauge theory approach. However, the existence 
of matter field can change the story completely. There 
is considerable experience in the study of a gauge field 
coupled to relativistic fermions and bosons, but little is 
known in the non-relativistic case. In the past few years 
there has been significant progress in the case of half- 
filling where there are no bosons. Consider the saddle 
point with ir flux through a plaquette. The fermionic 
spectrum are massless Dirac fermions with nodes at 
(jkj-,±Y-). Due to the additional SU(2) symmetry, 
these fermions are minimaly coupled to a set of SU(2) 
gauge fields. If the coupling to the gauge field is confin- 
ing, this leads to what is known in QCD as mass gener- 
ation and chiral symmetry breaking, which translates to 
gap formation and AF order in our language. [86] Thus 



the 7r flux phase is our route to AF order. On the other 
hand, consider the staggered flux state. This mean field 
ansatz breaks the SU{2) symmetry and the gauge field 
is broken down from SU(2) to U(l). This motivates the 
study of th problem of N 2-component massless Dirac 
fermions coupled to a U (1) gauge field in 2+1 dimensions. 
This model is often called QED3 in the literature. Since 
the staggered flux state has two spins and two nodes, the 
physical problem corresponds to N — 4. (It was shown 
that the velocity anisotropy of the Dirac cone is an irrel- 
evant variable and the low energy physics scales to the 
isotropic fixed or the QED3 model. [87|,[88j|) Assuming de- 
confinement, this class of states has been studied by 
expansion and is called the algebraic spin liquid. Her- 
mele et aZ.[90| showed, using results borrowed from the 
field theory literature, that for sufficiently large N, the 
deconfined state can be stable. In other words, the large 
N fixed point has no relevant perturbation, including the 
appearance of monopoles. This is at least a proof of the 
principle that deconfinement is possible in the presence 
of matter field. It is not known whether the critical N 
is greater than or smaller than 4, and here is where the 
QCD Monte Carlo community can help. (9^ 

Recently we received some encouraging news from ex- 
periments. It has long been thought that the Kagome 
lattice with S = \ has sufficient frustration to sup- 
port a spin liquid ground state. This seems to be re- 
alized by recent experiments on ZuCu3(OH)gCl2 where 
the Cu ions occupy Kagome sites and the system shows 
no AF order down to 30 mK, despite an AF exchange 
of ~ 200 K.[92j On the theoretical side, it is found that 
the projected wavefunction of a certain flux state gives 
the ground state energy obtained by exact diagonaliza- 
tion to within error. [93j This is remarkable for a trial 
wavefunction with no variation parameter, and is much 
better than what projected flux states did for the square 
lattice. The low energy physics of this state is that of two 
2-component massless Dirac fermions with spins, coupled 
to a U(l) gauge field. The confirmation of this picture 
for the Kagome lattice will give us reason to believe that 
the critical N is less than 4 in this case and perhaps in 
the case of the algebraic spin liquid as well. 

Prior to the Kagome example, there was strong evi- 
dence that a spin liquid state exists in the organic com- 
pound k-(ET) 2 Cu2(CN) 3 .[94|] Here the active sites are or- 
ganic molecular dimer molecules which form a triangular 
lattice to a good approximation. As mentioned in section 
III, the Heisenberg Hamiltonian on a triangular lattice is 
expected to order. Nevertheless, experimentally no AF 
is found down to 32 mK. The explanation lies in the fact 
that this system sits just on the insulating side of the 
Mott transition, so that charge fluctuations and ring ex- 
change terms lead to a more complicated effective Hamil- 
tonian which may favor the spin liquid state. [H, S Nu- 
merical work has given an adequate account of the phase 
diagram, including the appearance of superconductivity 
under pressure as shown in Fig. 3. [9^. l9a| On the trian- 
gular lattice the low energy Lagrangian is expected to be 
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an almost circular Fermi sea of spinons (fermions which 
carry S = i but no charge) coupled to U(l) gauge fields. 
This model has even more low lying fermionic excitations 
than the Dirac sea, and there is reason to believe that de- 
confinement occurs. Experimentally the observations of 
a constant spin susceptibility and linear specific heat at 
low temperatures are highly unusual for an insulator, and 
support the notion of a spinon Fermi surface. 



The appearance of particles, such as spinons which 
carry different quantum numbers compared with the orig- 
inal electron, is a phenomenon called "fractionalization." 
It is remarkable that the fermions originally introduced 
as a formal device in Eq. (4) take on a life of their own 
at the end of the day. Many people find this a hard con- 
cept to swallow and here is where explicit examples of 
a spin liquid will be a great help in convincing skeptics. 
I must emphasize that the new structure in the theory 
(spinons, gauge fields, etc.) emerges in the low energy 
physics and what emerges is independent of the way the 
problem was formulated. The constraint of no double 
occupation could have been enforced using Z(2), U(l) or 
SU(2) gauge fields, but the low energy structure will be 
the same. One particular formulation may simply be the 
convenient way to expose this emergent structure. 



The discovery of experimental examples of spin liq- 
uids is a very important development because for the 
first time, the deconfined (fractionalized) states are low 
temperature states which in principle can be studied in 
great detail. The success of the gauge theory method in 
these materials will give us confidence in its application 
to the more complex problem with holes. It is fortu- 
nate that the two promising examples of spin liquids are 
closely related to those discussed in the high T c context 
for underdoped and optimal doping, respectively. 



Finally, I would like to mention an approach to the 
pseudogap problem from a more general perspective. 
Senthil and I [§!| proposed that the underdoped cuprate 
should be considered as proximate to a spin liquid state 
which is then doped. In this view the pseudogap is the 
finite temperature region controlled by the quantum crit- 
ical point of doping a spin liquid by varying the chemical 
potential. The most promising spin liquid is in fact the 
algebraic spin liquid with its low energy massless Dirac 
fermion and U(l) gauge field. This point of view does 
not help solve the doped spin liquid problem, but does 
raise and answer the question: what is the signature of a 
deconfined state. The answer is that in the presence of a 
matter field, the only signature left is the irrelevance of 
instantons, i.e. the total gauge flux is conserved, much 
like the conservation of magnetic flux in our world. The 
experimental implication is subtle, but can be probed at 
least in principle in a suggested experiment. [99j] 



VII. QUANTUM OSCILLATIONS IN HIGH 
MAGNETIC FIELD 



Is there a region of the phase diagram where the pic- 
ture is less complicated and precise predictions can be 
made and tested? One natural idea is to apply a large 
magnetic field perpendicular to the layer to kill the su- 
perconductivity and ask what is the next most stable 
ground state which emerges. When H reaches H C 2 the 
vortex core overlaps and the core state becomes the ho- 
mogeneous ground state. Recent experiments indicate 
that we may b e probing this regime and have caused a lot 
of excit ement . 1 1 OOl 1 1 ll | It is believed that H C 2 for under- 
doped cuprates is approximately 100 T, well beyond what 
is available in the Laboratory. The resistive transition, 
on the other hand, is controlled by flux flow and is much 
more accessible to experiment. Recently, Shubnikov-de 
Haas oscillations are reported in an underdoped YBCO 
compound (T c = 55 K) in the flux flow region, in a mag- 
netic field range of 40 T to 62 T. From the period of os- 
cillation of the resistivity and Hall resistivity vs. -g the 
area of the Fermi surface pocket is extracted and found to 
contain 0.038 holes per pocket. If we assume the doubled 
unit cell scenario, there are two pockets in the reduced 
BZ and this corresponds to x = 0.076. The effective mass 
can also be measured from the temperature dependence, 
and m* = 1.9 m e . The doping concentration of YBCO 
is difficult to estimate, because the doping is from the 
oxygen chains and its charge state is not independently 
measurable. The best estimate by the experimentalists 
places x at 10%. More recently, similar quantum oscil- 
lations are reported in YBa2Cu40s, the so-called Y124 
compound (T c ~ 80 K) with magnetic field from 55 T 
to 85 T 10 1| and Shubnikov-de-Haas oscillations were 
reported in the range 45 T to 61 T. 102] The hole con- 
centration per pocket is now 0.05 and the inferred x of 
0.1 is again lower than the value of 0.125 that is believed 
to characterize this material. The effective mass is now 
3 m e . 



m 

Quantum oscillation is considered the best way to mea- 
sure the Fermi surface area, provided the sample is suf- 
ficiently free of defects that the electron can complete a 
cyclotron orbit without being scattered. Even at 50 T, 
the cyclotron orbit circumference is more than a thou- 
sand A in size (see Fig. 12), thus making it clear that 
disorder is really not an issue in the YBCO family. 

I think there is an excellent chance that the experi- 
ment is accessing the normal state that lies beyond H C 2- 
It is well known experimentally that quantum oscillations 
persist to H < H C 2- 103] The frequency of the oscillation 
is not changed upon crossing H c2 , only its magnitude 
is diminished. This phenomenon is particularly striking 
in layered materials where the flux flow regime can be 
very wide. For example, a recent experiment on organics 
observed quantum oscillations down to \H C 2 and is even 
more striking in that the exp ecte d reduction of the ampli- 
tude was not even observed. jl04j A rough physical expla- 
nation may be as follows. Let us take a cartoon picture 
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FIG. 12: Vortex cores are represented by discs of radius R v . 
Figure is drawn almost to scale at H = \H C 2, when the discs 
occupy half of the total area (R v « 25 A). A cyclotron orbit 
is shown, also roughly to scale. Its area encloses about 10 
vortices for the experiment on YBCO. The long dimension of 
the orbit is about 800A. 



of the vortex core as discs of radius R v . At H C 2 the discs 
overlap and we define H C 2 = 4>q/(ttR 2 ) where 4>o = hc/2e 
is the superconducting flux quantum. Note R v = \/2£ 
where £ is the coherence length. For H C 2 — 100 T, R v 
— 25 A. For H = \H c2 the disc density is rather high, 
as drawn in Fig. 12, almost to scale. Let us estimate the 
size of the semiclassical orbit. Very generally, the real 
space area A is related to the momentum space pocket 
size Ak by 



latter case, the I vector in all the vortices inside the cy- 
clotron orbit are either pointing up or pointing down, i.e. 
the coherence length of the order must exceed the orbit 
size. For the aspect ratio chosen, the long dimension of 
the orbit is about 800 A. What is required for quantum 
oscillations is strong tunnelling between the vortex cores. 
Then the quasiparticles are extended states which carry 
information of the uniform normal state for H > H C 2 be- 
cause the random phase coming from tunnelling through 
the superconducting region will cancel. Given the dense 
packing of vortices, this seems quite reasonable. 

What is the origin of the pockets that give rise to the 
quantum oscillations? Here we are on rather uncertain 
grounds. In general, we can classify the state as conven- 
tional or exotic. By conventional I mean there is some 
new order giving rise to unit cell doubling and a Fermi 
liquid ground state, and exotic means ever ythin g else. 
While exotic scenarios have been suggested [105} here I 
will concentrate on discussing conventional possibilities, 
knowing that we immediately encounter a problem con- 
cerning the relationship of the measured area to doping 
estimates. Either the doping estimate is incorrect, or 
there exists undetected electron or hole pockets. Let me 
forge ahead and assume that the doping concentrations 
are actually x = 0.075 and 0.10 and correspond to the 
pocket areas measured in the two experiments and dis- 
cuss several candidates for the unit cell doubling order. 



A = A k 



he 



(15) 



It is useful to introduce the Landau length 



f-H = \ — 77- 



(16) 



Note that the flux through an area n£ 2 H is </>o, so that 
R v = £fj=H a2 ■ The Landau level number v at the Fermi 
level is given by the number of full flux quantum cf>o = 2</>o 
which penetrates A, i.e. 
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We find v = 10 for H = 50 T, using the experimentally 
measured A k = 5.1 x 10 14 cm~ 2 for YBCO. Since 50 T 
corresponds to « ^H c2 , this also means that on average, 
10 vortex cores are inside the spatial orbit. Since the cy- 
clotron motion is fast compared with that of the vortices, 
we can assume a static picture of the vortex core, which 
either forms a hexagonal lattice or somewhat distorted 
from it. Assuming the orbit to be an ellipse with aspect 
ratio ~ 8 (for reasons to be explained later), we sketch 
a snapshot of the orbit in relation to the vortex core 
in Fig. 12. In order to see quantum oscillations, what 
is required is that the states in the cores are in phase, 
whether it is AF order or staggered-flux order. In the 



AF order. We saw in section I that at very low dop- 
ing, hole pockets appear around (n, n) in the back- 
ground of AF order. The question is whether this 
narrow region (x < 3%) will open up in a strong 
magnetic field and extend up to 10 or 12% doping. 
This is probably what most people have in mind 
when they refer to AF in the core and we shall call 
states of this kind "doped AF." I think this sce- 
nario is unlikely. From neutron scattering we know 
that there is a sharp triplet resonance at ~ 30 meV 
in this doping range at (i t, tt) and a sharp drop 
in spectral weight below it.[107[ AF order requires 
the condensation of this triplet mode. A 50 T field 
corresponds to an energy scale of 6 meV, hardly 
enough to perturb the resonance mode. Thus I 
think the AF is energetically highly unfavorable. 
Nevertheless, there are experimental ways to dis- 
tinguish doped AF from the staggered flux state, 
which I will return to later. 

Staggered flux state. As we argued before, this 
state is energetically favored to appear once super- 
conductivity is destroyed. I want to point out that 
in principle a small amount of AF order can co-exist 
with this state. As seen in Fig. 8, this is favorable 
for x < 0.08. However the driving force and excita- 
tion spectrum of this state is completely different 
from the doped AF. In that case the holes live in 
the lower Hubbard band, i.e. the band is a down- 
ward parabola, separated by a gap of order of the 
Mott gap with energy scale 1 eV. 
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Valence bond solid (VBS), or more general states 
which break translation symmetry without break- 
ing time reversal as in AF. There has been dis- 
cussion of this state living in the vortex core. jl06j 
Roughly speaking, the electronic spectrum may be 
quite similar to the staggered flux state except that 
it is expected to be gapped. 



Spin density wave (SDW). [1081 ] This seems very 
attractive phenomenologically. From NMR there 
are c laim s for enhanced antifcrromagnetism in the 
core. [l0!| llldf l In the doped La2CuC>4 system it 
is known from neutron scattering that incommen- 
surate S DW is in duced in the vicinity of the vor- 
tex core. [lllL Ill2| This is quite understandable be- 
cause inelastic neutron scattering indicates soft ex- 
citations at precisely these incommensurate wave 
vectors, typically around |E from (f,f)- This 
is often connected to the notion of stripes near 
the concentration x = |. However, static SDW 
has never been seen in the YBCO family. Neu- 
tron scattering has reported dynamical fluctuations 
at 20 meV at the incommensurate wave vectors 
(t?- ± S, 7?- ± S) . While recent data on more under- 
doped samples (T c = 50 K) sees dyamic scattering 
down to 12 meV, it has much less spect ral weight 
than in the doped La 2 Cu04 system. [l07j Thus it is 
questionable whether a magnetic field can induce 
SDW order. Furthermore, the prevailing view is 
that the SDW is not tied to Fermi surface nesting, 
and it is unclear whether hole pockets should be 
expected. 

How can we distinguish the staggered flux state 
from the doped AF state? A useful experimental 
test is to do a tunnelling experiment to measure the 
density of states in high fields. Since spatial reso- 
lution is not needed, we do not need STM, which is 
probably impossible under these pulsed field con- 
ditions, but we propose either break-junction tun- 
nelling or other tunnelling geometry where elec- 
trons can be injected into the a-b plane from the 
side. Schematically the density of states of the 
staggered flux state vs. the doped AF is sketched 
in Fig. 13(a). Of course, we expect much of the 
higher energy features to be smeared out, but 
strong particle-hole asymmetry and a dip in the 
density of states above Sf will support the stag- 
gered flux picture. On the other hand, if this state 
is the continuation of the doped AF state observed 
at x < 0.03 for zero field, we expect a large Mott 
gap above the Fermi energy, which has been esti- 
mated to be 30 meV in the case of Y124.[l0l| Thus 
tunnelling data is a sensitive test of the notion of 
hole pockets in the lower Hubbard band. 

Next we extract some physical parameters from the 
data. Let us parametrize the Dirac spectrum of the 
staggered flux state near each node by 
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FIG. 13: Schematic picture of the tunneling density of states 
for (a) the spin split bands in the doped staggered flux state 
and (b) the doped AF. If the hole pocket lies in the lower 
Hubbard model, a large Mot gap is expected. 
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where k\{k2) in the direction perpendicular (paral- 
lel) to the large Fermi surface. By definition 
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By writing the density of states m* /2n in terms 
of the velocities along the Fermi surface, we find 
the anisotropic generalization of the relation vp — 
htcF/m : 



(20) 



Equation (20) applies equally well to the 
anisotropic Dirac spectrum and the anisotropic 
parabolic band. The right hand side of this equa- 
tion are directly measured. Putting in the numbers 
for YBCO with m* = 1.9 m e , we find hy/vpuK 
to be 0.48 eVA. The Fermi velocity can be mea- 
sured by ARPES, even though the error bar is sub- 
stantial. Taking the value vf = 1-4 eVA from 
Campuzano et aL,[3l[ we find the anisotropy to 
be vf/va — 8.4. This compares with the value 
7.9 directly extracted from thermal conductivity 
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data ^8l|] in their T c — 62 sample. In that pa- 
per they extract a gap value A of 71 meV, using 
vp = 1.65 eVA. With our numbers we expect Ao = 
57 meV, in very good agreement with ARPES esti- 
mate of the pseudogap for a similarly doped Bi 2212 
sample. A similar exercise using Y124 data [lOl| 
yields vf/va — 16.7 and A = 29 meV. The gap 
value is on the low side. Given that these esti- 
mates are sensitive to the square of m* and Vf and 
their associated errors, these are reasonable num- 
bers. More importantly, it shows the correct trend 
of increasing anisotropy and reduced gap with in- 
creasing doping. It is worth noting that the el lipse 
drawn in Fig. 1(b) by Doiron-Leyraud et aZ. |l00l | 
has an aspect ratio of vf/va ~ 4. Noting that the 

length of the ellipse scale as y~%^> our numbers 

indicate a more elongated ellipse and a small de- 
parture of the inside edge of the pocket away from 
the Fermi arc. For Y124, the ellipse almost reaches 
the saddle point at (0, 7r). Beyond this point the 
pockets connect to form large Fermi surfaces which 
are presumably unobservable due to disorder scat- 
tering, and we expect the quantum oscillation to 
disappear. There will be a transition to the uni- 
form large Fermi surface (area 1 — x) state, but 
the exact nature of this transition is an open issue. 
We reiterate that Eq. (20) applies to the parabolic 
band of the doped AF as well. The only difference 
is that in that case there is no reason to identify vp 
with the normal state Fermi velocity. 

Quite generally, if this state is a Fermi liquid in a 
nonmagnetic background, we expect Zeeman split- 
ting of the up- and down-spin pockets. For free 
fermions the splitting is entirely determined by m* 
which gives the density of states. Let us define 
AAk = Afef — Ak[ to the difference in the up- and 
down-spin pockets. We find 



AA k = A£— = AE2Trm*/h 2 (21) 
dE 

Using AE = g[i B H with g = 2 and H = 60 T, 
we find for Y124 a splitting of AA k /A k m 26%. 
Experimentally a splitting of w 10% was observed. 
If we attribute this to spin splitting, we will need 
a reduction of the predicted magnetization due to 
residual interactions. In Fermi liquid theory, this 
is given by the Landau parameter Ffi, so that x — 
Xo/(1 + Fq). Thus we need Fg = 1.6 which is 
quite reasonable for a strongly interacting system. 

In contrast, for the doped AF, we do not expect 
a similar spin splitting. This was pointed out to 
me by T. Senthil. This is because the AF will cant 
with sublattice magnetization perpendicular to the 
applied field and spin along the field is no longer 
a good quantum number. If splitting due to inter- 
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FIG. 14: Schematic phase diagram predicted for an under- 
doped cuprate in a magnetic field. The solid lines are true 
phase transitions while the dashed lines are cross-overs to the 
pseudogap phase and the vortex liquid phase. The vortex 
liquid at H = corresponds to the Nernst region in Fig. 2. 



layer tunnelling can be ruled out, this is another 
argument against the doped AF. 

Finally, we discuss the cyclotron resonance exper- 
iment. In principle, this can help distinguish be- 
tween the Dirac spectrum of the staggered flux 
state and the parabolic spectrum expected for the 
doped AF.(7f| As made fam ous by recent experi- 
ments on graphene. [ll3l [rTij the Landau level in a 
Dirac spectrum is given by E v = \ 2vH-JvfvaI ^-bh 
which scales as \[B, instead of the standard Landau 
level spacing heB/m*c which is linear in B in con- 
ventional parabolic bands. Using measured param- 
eters for YBCO we find the splitting between the 
v = 9 and v = 10 Landau level to be 3.1 meV, com- 
pared with 3.3 meV for a parabolic band, an un- 
fortunate coincidence. In principle, one can check 
the \/B vs. B dependence, but this may be a chal- 
lenging experiment to carry out in pulsed fields. 

To summarize, with insight gained from Nernst ef- 
fect measurements, we propose in Fig. 14 a phase 
diagram for the underdoped cuprate in a field. The 
staggered flux state appears as a new phase in high 
field and low temperature. Since this state requires 
the phase coherence of the bosons just as the su- 
perconducting state, we expect its transition tem- 
perature to be comparable to the superconducting 
T c . We expect a signature in transport properties 
at this transition temperature. The H C 2 line is a 
cross-over line which characterizes the onset of a 
vortex liquid state. This is the extension of the 
Nernst region shown in Fig. 2 to finite H. The true 
superconducting transition is to a vortex solid or 
glass state. Whether the staggered flux state ter- 
minates in the vortex liquid or solid phase is not 
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known. Below a large energy scale Ao lies the spin 
gap phase which is the extension of the pseudogap 
phase in Fig. 2 to finite H . We emphasize that 
while the model of Hsu et and the more phe- 

nomenological approach of Chakravarty et aL[69| 
also predict a staggered flux (DDW) state at high 
field, they key point they made is that this state 
is smoothly connected to the pseudogap phase, i.e. 
the solid phase boundary in Fig. 14 is pushed up to 
an energy scale A . Great effort has been made to 
look for the moment generated by the orbital cur- 
rents by neutron scattering and nothing convincing 
has been found. There is no sign of a phase tran- 
sition even in the highest purity materials. It is 
difficult to see how the Fermi arc can be consistent 
with a Fermi liquid picture of sharp quasiparticles. 
In my opinion, existing data at low field is suffi- 
cient to rule out the DDW scenario. In contrast, in 
our picture the pseudogap phase is separated from 
the high field staggered flux Fermi liquid state by 
a phase boundary. While a detailed description of 
the Fermi arc is not available, it is at least possible 
to imagine that different physics is at play and the 
Fermi arc picture can be consistent with the exper- 
imental observation of quantum oscillations in high 
fields. 

Finally, the proposed phase diagram is heavily 
based on mean field theory and the notion of bo- 
son condensation and should be viewed simply as 
a starting point for discussion of the experiment. 
It is possible that the orbital current can survive 
fluctuations, but it is also possible that tunnelling 
events wipe out the ordering pattern, leaving us 
with an exotic ground state. Other order, such as 
weak AF, could also emerge to take its place. The 
overall picture, I think, is more general than the 
survival of the orbital currents. The large spin gap 
at (0, 7r) and features in the tunnelling density of 
states sketched in Fig. 13(a) may survive. Clearly 
more experiments in high fields will be needed to 
shed light on these important issues. 



VIII. QUESTIONS AND ANSWERS 

Instead of a conclusion, I would like to pose a number 
of questions and attempt to answer them. 

1. What is unique about the cuprates? What is the 
pairing mechanism? What controls T c ? Can T c be 
raised? 

The cuprate is the prime example of a two- 
dimensional Mott insulator which can be doped. 
The copper 2+ ion is in a d 9 configuration, making 
it a spin i system. The layer symmetry breaks the 
degeneracy of the d orbitals dow to a single d x 2_ y 2 
orbital, so that orbital issues which cause complica- 
tions in most other transition metal oxides do not 



play a role here. In most other examples, strong 
coupling of the charge state to distortion of the 
oxygen cage causes doped carriers to localize. For 
example, if we replace Cu by Ni, La2Ni04 requires 
50% hole doping before it becomes metallic. In 
many ways, the cuprate is the ideal model system 
for a doped Hubbard model. The parent compoud 
is one of the few layered S = ^ antiferromagnet 
with no orbital degeneracy that is known. It is 
hard to find a material that is so simple! 

Searching for a pairing mechanism may be the 
wrong question, because implicitly the assump- 
tion is made that pairing occurs via the conven- 
tional route of well defined quasiparticles exchang- 
ing some bosonic glue. 115| We know that under- 
doped cuprates do not fit this description. In- 
stead superconductivity is the ground state which 
emerges as the best compromise between the ki- 
netic energy of hole motion and the antiferromag- 
netic exchange. The energy gap is set by the ex- 
change energy J while T c is controlled by phase co- 
herence and is of order xt. Thus the optimal T c is 
set by J. The cuprate family already has the high- 
est J among transition metals. (The ladder com- 
pound made up of the same copper-oxygen bond is 
the only example I know which has a higher J than 
the planar cuprate.) The hopping matrix element 
t is also large because of the strong covalent bond- 
ing between the d x 2_ y 2 orbital with the oxygen p 
orbital. Thus the cuprates are blessed with a large 
energy scale. The variation of T c on a few tens 
of degrees between YBCO and the Hg compound 
seems to be attributable to changes in the further 
neighbor hopping matrix elements. [l9| 

If the superconductivity is due to the electronic en- 
ergy scale such as J and t, why isn't T c higher? 
As we saw in section I, it is only ~ ^ or ^ for 
both cuprates and organics. What accounts for 
the small numerical factor? One consequence of 
doping a Mott insulator is that the carrier density 
and therefore the superfluid density is small. Thus 
phase fluctuation is important and limits T c . In a 
d-w&ve superconductor, the thermal excitation of 
nodal quasiparticles further suppresses the super- 
fluid stiffne ss on a short distance scale, as seen by 
experiment. [HH This reduces the T c further. This 
is why T c is such a small fraction of J. Besides 
increasing the overall energy scale, one possible av- 
enue for enhancing T c is to try to stabilize super- 
conductors with a full gap such as d+id states. For 
example, on a triangular lattice, the d + id state is 
believed to be favored. In this case, the quasipar- 
ticle suppression of the superfluid density is negli- 
gible and T c will be higher, everything else being 
equal. 

2. What about other approaches? Why is there no 
consensus after all these years? 
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In a review of this kind which adopts a particular 
viewpoint (going under the names strong correla- 
tion, RVB, gauge theory, etc.) I feel obligated to 
give my honest opinion on other approaches. This 
may not be the best way to win friends, but let 
me proceed. First let me point out that within the 
RVB approach there are many variants that have 
not been covered here. I emphasized the formu- 
lation based on a fermionic representation of the 
spins and a bosonic representation of the holes. The 
opposite formulation with bosonic spinons [l06j | 
and a mo re co mplicated formulation with flux at- 
tachment [l 1 7l ] has also been pursued. Below let me 
focus on approaches outside the RVB "big tent" by 
separating the most prominent work into four ma- 
jor categories and discussing them in turn. 

(a) Spin fluctuation models. 

Historically the possibility of <i-wave super- 
conductivity based on the exchange of AF spin 
fluctuations was discussed prior to the discov- 
ery of high T c and a numb er of workers apply 
these ideas to the cuprates. [TTj. Ul l Il20l . fl2l| 
These discussions are either based on random 
phase approximation treatment for the Hub- 
bard model, or phenomenological coupling be- 
tween spins and electrons treated as sepa- 
rate entities. The latter may be correct in 
heavy fermion systems, where there indeed ex- 
ists separate itinerant and local moments, but 
is more questionable in the Hubbard model 
where there is only one kind of electronic state. 
The best chance for this model to work is 
in the overdoped region, where well defined 
quasiparticles exist in the normal state and 
the onset of superconductivity appears quite 
conventional. However, from neutron scatter- 
ing it is known that spin fluctuations become 
very weak, with spin correlation length of 10 A 
or less beyond optimal doping. A variant of 
this idea ascribes the boson being exchanged 
to the sharp triplet resonance discovered by 
neutron scattering. However, this resonance is 
very narrow in momentum space and carries 
very little spectral weight, ma king it unlikely 
to be the powerful glue needed. [l22j | It is clear 
that in order for the spin fluctuation picture 
to work, one has to deal with the strong cou- 
pling problem and include spin fluctuations at 
all energy scales up to J. Numerical meth- 
ods may be the only option. Recent numerical 
studies based on cluster dynamical mean field 
theory shows promise in this direction in that 
a pairing kernel was identified in the triplet 
channel which peaks at and resembles 

the spin fluctuation channel. [53| I think the 
most serious limitation of this line of ideas is 
that it completely fails to address the pseudo- 



gap issue. As we have seen, there is simply 
no well defined quasiparticle above T c to pair 
them by whatever mechanism. By avoiding 
the issue of proximity to a Mott insulator, this 
approach misses the key physics and the most 
interesting aspect of high T c , 

(b) Microscopic inhomogeneity, stripes, etc. 

The basic idea is that one way the doped Mott 
insulator can resolve the competition between 
kinetic energy and AF exchange is for the 
holes to phase segregate. In particular, the 
holes may concentrate into one dimensional 
regions where AF order is suppressed, sepa- 
rated by an AF regio n where the hole den- 
sity is suppressed. [123| Experimentally, it was 
discovered that a particularly stable configu- 
ration where the holes form a quarter filled 
chain (one hole per two sites along the chain) 
separating AF regions as an anti-phase bound- 
ary ind eed exists in L^-^Ba^CuC^ where 
x — 7T. [l24| This configuration and variants 
thereof has been called stripes. The stripe 
explains the suppression of T c in the doped 
La2CuC>4 system near x = It also re- 
ceives theoretical support from numerical so- 
lution of the t-J model using density ma- 
trix renormalization group method. ,52] Thus 
there is good reason to believe that stripes 
form a strongly competitive ground state, es- 
pecially near x = g. However, the order- 
ing temperature for the charge order which 
preceded spin order is low, of order 40 K, 
and the spin order (incommensurate SDW) 
is even lower. Furthermore, static stripes are 
seen only in the doped La2Cu04 family and 
not in the YBCO family, and in the latter 
case, even dynamical stripe fluctuations have 
weak intensities, as discussed earlier. Thus 
it is my belief that stripes represent low en- 
ergy physics and emerge as the ground state 
only as a result of delicate competition. Re- 
cently there is a remarkable report that the 
ARPES in Lai. 8 7 5 Ba.i25Cu04 reveals quasi- 
particles which d-wave dispersion with a max- 
imum gap of 20 meV, as if the existence of 
stripes has no effect on pairing of the pseu- 
dogaps except at very low energy. [125J This 
reinforces my view that stripe order is a sec- 
ondary effect. 

Static stripes are clearly detrimental to super- 
conductivity. Nevertheless, there has been a 
strong push, notably by S. Kivelson and co- 
workers, for the idea that fluctuating stripes 
may be responsible for superconductivity and 
the pseudogap phenomenon. I have always 
found the motivation of this idea puzzling. 
Fortuna tely a recent article by Kivelson and 
Fradkin jl26l | lays out the physical ideas and 
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motivation in a very clear way. They share our 
basic assumption that doping a Mott insulator 
is the key question and they share our notion 
that superconductivity comes from strong re- 
pulsion. However, they reject the possibility 
that this can happen in a homogeneous phase. 
As far as I can tell, the reason offered is that 
"the enormous effort has been devoted to nu- 
merical searches for superconductivity in vari- 
ous uniform Hubbard and t-J related models, 
with results that are, at least, ambiguous" and 
the feeling that if that were the right route, 
"unambiguous evidence of it would have been 
found by now." Instead they focused on mod- 
els where superconductivity from a repulsive 
Hubbard model can be demonstrated, namely 
the ladder system. The ladder system can be 
considered an example of a valence band solid, 
where spins on the rung form singlets and 
doped holes tend to form pairs on the rung. 
Then pair hopping between ladders provides 
bulk superconductivity. This is all well and 
good but Kivelson and Fradkin are the first 
to admit (in fact insist) that this toy model 
has little to do with cuprates. For example, 
it is hard to imagine how quasiparticles can 
traverse this system or ladder (or any kind of 
stripes) at a 45° angle and become the gap- 
less nodal particle. As long as this problem 
is not addressed adequately, I fail to see how 
this approach can be a viable route to super- 
conductivity in the cuprates. 

(c) Phonons. 

Electron phonon coupling is typically strong 
in transition metal oxides involving the E g or- 
bitals, because of strong coupling to the dis- 
tortion of the oxygen cage. As mentioned ear- 
lier, in most cases this leads to localization of 
the doped hole. While holes in the cuprate es- 
cape this fate, the electron phonon coupling 
must be strong. Nevertheless, we have ar- 
gued that the dominant energy in this problem 
is Coulomb repulsion. Hence the discussion 
of electron-phonon coupling must be made in 
this context. A combination of numerical and 
analytic work has shown that short range re- 
pulsion has the effect of reducing electron- 
phonon scattering at large momentum trans- 
fer because charge fluctuati ons at short dis- 
tance is suppressed. [l27l[l28j It is also claimed 
that under certain circumstances, d-wave pair- 
ing is possible with electron-phonon coupling, 
contrary to conventional wisdom. Experimen- 
tally, phonon sidebands appear visible in STM 
tunnelling spectra. All this says that electron- 
phonon coupling is there (no great surprise) 
but does not provide any evidence that it is 
the essential driving force behind pairing. 



The isotope effect is often quoted as evidence 
for the phonon mechanism for high T c . It 
turns out that isotope effects affect T c via 
the superfhhd density, and not the order pa- 
rameter itself. [129l | Thus the isotope effect in 
fact supports the view that T c is controlled by 
phase fluctuations in the underdoped region. 
There is no isotope effect at optimal doping 
and beyond. 

(d) Three band model. 

Since the copper oxygen layer involves one 
copper and two oxygen per unit cell (excluding 
the apical oxygen), the minimal microscopic 
electronic model requires a d x 2_ y 2 Cu orbital 
and two oxygen p or bita ls. This is called 
the three band model. [l30l Il3l| Various clus- 
ter calculations indicate that the low energy 
physics (below a scale of f ~ 4 meV) can be 
adequately described by a one band Hubbard 
model. [H Over the years this has become the 
majority view, but there are still workers who 
believe that the three band model is required 
for superconductivity. I think the original mo- 
tivation for this view came from a period of 
nearly ten years when much of the community 
believed that the pairing symmetry is s-wave. 
It certainly is true that a repulsive Hubbard 
model cannot have an s-wave superconducting 
ground state. Three band models with further 
neighbor repulsion were introduced to gener- 
ate the requisite effective attraction. The s- 
wave story was overturned once and for all 
by phase sensitive experiments in 1994 but 
some of the three band model proponents per- 
severed. In particular, Varma introduced the 
idea of intra-cell orbital currents, i.e. a cur- 
rent pattern flowing between the Cu-0 and 
0-0 bonds as a model for the pseudogap.[132j 
This has the virtue of leaving the unit cell in- 
tact and this kind of order is very difficult to 
detect. With such a complicated model it is 
difficult to make a convincing case based on 
theory alone and a lot of attention has been 
focused on experimental detection of time re- 
versal symmetry breaking or spontaneous mo- 
ments due to orbital currents. Unfortunately, 
the expected signals are very small and can 
easily be contaminated by a small amount of 
minority phase with AF order. Up to now 
there is no reliable results in support of this 
kind of orbital currents. 

Why is the high T c problem hard? Why is the high 
T c problem important? 

One of the ironies about the high T c problem is that 
the ground states are rather conventional. It ranges 
from AF to d-wave superconductors to Fermi liq- 
uid as the doping is increased. The exciting new 
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physics happens at finite temperatures where Fermi 
arcs, pseudogaps and other novel phenomena ap- 
pear. This region is not characterized by order in 
the traditional Landau sense, and it is not possible 
to make precise statements as one could if one has 
an exotic ground state. From the theoretical side, 
what we can do so far to classify the pseudogap as 
the deconfined phase of a gauge theory, where new 
particles, spinons and holons and £7(1) gauge fields 
emerge in the low energy physics. These particles 
and gauge fields are still strongly coupled and we 
are not able to compute and compare with experi- 
ments in a quantitative way. We can make various 
caricatures (such as local boson condensation and 
fluctuations between different phases) which mimic 
the observations in a qualitative manner. Neverthe- 
less, I believe the problem is so strongly constrained 
by data that this path is the only viable option at 
the moment. I must emphasize that the gauge field 
is an integral part of this story, and not a technical 
artifact of our particular formulation. In principle 
it is possible to formulate the problem with spin 
carrying bosons and fermions and holes, and reach 
the same low energy state described here. In this 
sense this proposed solution of the high T c problem 
represents the emergence of new physics. 

Two recent developments are encouraging. First, 
it appears possible to probe the "normal state" 
at zero temperature by applying a large magnetic 



field. While the state again appears to be a Fermi 
liquid, it may expose a new kind of order and re- 
veal what is fluctuating in the pseudogap phase. 
Secondly, the elusive spin liquid ground state may 
finally have been realized. Even though this is an 
early stage in its development, the study of the low 
energy excitations of these ground states may fi- 
nally allow us to access spinons and gauge fields. 
The effective low energy theory of the spin liquid 
problem is also more amenable to treatment be- 
cause bosonic holes are not involved. It is our view 
that the study of spin liquids and high T c opens a 
new chapter in condensed matter, one that defies 
the dual Landau paradigm (local order parameter 
and Fermi liquid theory) that has been the back- 
bone of our Science for three quarters of a century. 
To me, this is why the high T c problem is important 
and exciting. 
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